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Abstract 

We exhibit conditions under which the flow of marginal distributions 
of a discontinuous semimartingale ^ can be matched by a Markov process, 
whose infinitesimal generator is expressed in terms of the local character- 
istics of ^. Our construction applies to a large class of semimartingales, 
including smooth functions of a Markov process. We use this result to de- 
rive a partial integro-differential equation for the one-dimensional distri- 
butions of a semimartingale, extending the Kolmogorov forward equation 
to a non-Markovian setting. 
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1 Introduction 



Stochastic processes with path-dependent / non-Markovian dynamics used in 
various fields such as physics and mathematical finance present challenges for 
computation, simulation and estimation. In some applications where one is 
interested in the marginal distributions of such processes, such as option pricing 
or Monte Carlo simulation of densities, the complexity of the model can be 
greatly reduced by considering a low-dimensional Markovian model with the 
same marginal distributions. Given a process ^, a Markov process X is said to 
mimick ^ on the time interval [0, T], T > 0, if ^ and X have the same marginal 
distributions: 

Vfe[0,T], ^t^Xt. (1) 

X is called a Markovian projection of i^. The construction of Markovian projec- 
tions was first suggested by Bremaud ^ in the context of queues. Construction 
of mimicking processes of 'Markovian' type has been explored for Ito processes 
p2] and marked point processes f^. A notable application is the derivation of 
forward equations for option pricing [31 [3] . 

We proposer in this paper a systematic construction of such Markovian pro- 
jections for (possibly discontinuous) semimartingales. Given a semimartingale ^, 
we give conditions under which there exists a Markov process X whose marginal 
distributions are identical to those of ^, and give an explicit construction of 
the Markov process X as the solution of a martingale problem for an integro- 
differential operator [1 [H [H [H]. 

In the martingale case, the Markovian projection problem is related to the 
problem of constructing martingales with a given fiow of marginals, which dates 
back to Kellerer |18j and has been recently explored by Yor and coauthors 
[U [131 [2D] using a variety of techniques. The construction proposed in this pa- 
per is different from the does not rely on the martingale property of ^. We shall 
see nevertheless that our construction preserves the (local) martingale property. 
Also, whereas the approaches described in [H [M] [20] use as a starting point the 
marginal distributions of ^, our construction describes the mimicking Markov 
process X in terms of the local characteristics Il6| of the semimartingale 
Our construction thus applies more readily to solutions of stochastic differen- 
tial equations where the local characteristics are known but not the marginal 
distributions. 

Section[2]presents a Markovian projection result for a R'^-valued semimartin- 
gale given by its local characteristics. We use these results in section l^T^ to derive 
a partial integro-differential equation for the one-dimensional distributions of a 
discontinuous semimartingale, thus extending the Kolmogorov forward equation 
to a non-Markovian setting. Section [3] shows how this result may be applied 
to processes whose jumps are represented as the integral of a predictable jump 
amplitude with respect to a Poisson random measure, a representation often 
used in stochastic differential equations with jumps. In Section [4] we show that 
our construction applied to a large class of semimartingales, including smooth 
functions of a Markov process (Section [4. ip . and time-changed Levy processes 
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(Section 1121) . 



2 A mimicking theorem for discontinuous semi- 
martingales 

Consider, on a filtered probability space (fi, J", (J-t)t>o,P), an Ito semimartin- 
gale, on the time interval [0,T], T > 0, given by the decomposition 

^t^^o+ f Psds+ f SsdW,+ f f yM{dsdy)+ j f yM{dsdy), 

Jo Jo Jo J\\y\\<l Jo J\\y\\>l 

(2) 

where is in R'', is a standard M" -valued Wiener process, M is an integer- 
valued random measure on [0, T] x with compensator measure /i and M = 
M — /i is the compensated measure [TH Ch.II,Sec.l], (3 (resp. 5) is an adapted 
process with values in M'' (resp. Mdxn(R))- 

Our goal is to construct a Markov process, on some filtered probability space 
(J7o, B, {Bt)t>o, Q) such that X and ^ have the same marginal distributions on 
[0,T], i.e. the law of Xt under Q coincides with the law of under P. We 
will construct X as the solution to a martingale problem [TTJ [531 US HT] on the 
canonical space = D([0, T], R'*). 



2.1 Martingale problems for integro-differential operators 

Let — Z?([0, T], K'^) be the Skorokhod space of right-continuous functions 
with left limits. Denote by Xt{io) — Luit) the canonical process on rig, its 
filtration and Bt = B^_^. Our goal is to construct a probability measure Q 
on fio such that X is a Markov process under Q and ^ and X have the same 
one-dimensional distributions: 

Vte[0,T], it^Xt. 

In order to do this, we shall characterize Q as the solution of a martingale 
problem for an appropriately chosen integro-differential operator C. 

Let C^{W^) denote the set of bounded and continuous functions on R"^ and 
C^(R'') the set of infinitely differentiable functions with compact support on R''. 
Consider a time-dependent integro-differential operator C = {Ct)te[o,T] defined, 
for / G C^^(R''), by 

«,i=i ' ^ (3) 

+ / [f{x + y) - f{x) -l{i\y\\<i}y.Vf{x)]n{t,dy,x), 



where a : [0, T] x R'^ i-^ A/<jxd(R), & : [0, T] x R'^ i-^ R"* are measurable functions 
and, for each {n{t, . ,.t), {t,x) G [0,r] x R'') is a measurable family of positive 
measures on R'' — {0}. 
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For xq in R'^, we recall that a probability measure Qxo on {il,o,BT) is a 
solution to the martingale problem for (£, C5"(R'*)) on [0, T] if Q (Xq = xq) = 1 
and for any / e C(J°(R''), the process 



fiXt) - /(a:o) 



/ Csf{Xs)ds 
Jo 



is a (Qj^Q, (Z?()t>o)-martingale on [0,r]. Existence, uniqueness and regularity 
of solutions to martingale problems for integro-differential operators have been 
studied under various conditions on the coefficients [151 HH [121 HI] ■ 
We make the following assumptions on the coefficients: 

Assumption 1 (Boundedness of coefficients). There exists Ki > 

V(i,z) e [0,T] xM'^ \\bit,z)\\ + \\ait,z)\\+ f {l A \\yf) n{t,dy., z) < 



and 



lim / sup n{t,{\\y\\> R},z) dt = 0. 



where \\.\\ denotes the Euclidean norm. 

Assumption 2 (Continuity). For all t e [0,T] and B e B{R'^ - {0}), b{t,.), 
a{t, .) and , n{t, B, .) are continuous on W^, uniformly in t ^ [0, T]. 

Assumption 3 (Non-degeneracy). 

Either \JR > G [0, T] inf inf *a:.a(i, z).x > 

\\z\\<R £ceR<i 

or a = and there exists/3 s]0, 2[, C > 0, and a family n^{t, dy, z) 
of positive measures such that 

V(t, z) G [0, T] X M." nit, dy, z) = (i, dy, z) + dy, 



(lA||y||^ 



n^{t,dy,z) < K2, lim sup / Hyf n''(i, dy, z) = 0. 



Mikulevicius and Pragarauskas [2T] show that if C satisfies Assumptions 1-3 
(in which corresponds to a "non-degenerate Levy operator" in the terminology 
of [21]) the martingale problem for (£, Cg^{M.'^) ) has a unique solution for every 
initial condition xq G M'': 

Proposition 1. 121, Theorem 5] 

Under Assumptions[Il[Band\^ the martingale problem for {{Ct)te[Q.T],Ca^{^'^)) 
on [0, T] is well-posed : for any xq G , there exists a unique probability measure 
Q^o on {D[[Q,T],W^),Tt) such that Qx„{Xo = xo) = 1 and for any f e Cf^iR'^) 

fiXt)-fixo)- f Csf{Xs)ds 

JQ 



4 



is a Qa;p -martm^aZe. Under Q^o; i^t) is o, Markov process and the evolution 
operator (Qt)te[o,T] defined by 

V/eCO(R'^) Qt/(xo) = EQ-o [/(Xt)] (4) 
is strongly continuous on [0, T[. 

2.2 A uniqueness result for the Kolmogorov forward equa- 
tion 

An important property of continuous-time Markov processes is their link with 
partial (integro-) differential equation (PIDE) which allows to use analytical 
tools for studying their probabilistic properties. In particular the transition 
density of a Markov process solves the forward Kolmogorov equation (or Fokker- 
Planck equation) . The following result shows that under Assumptions [TJ [5] 
and |3] the forward equation corresponding to C has a unique solution: 

Theorem 1 (Kolmogorov Forward equation) . Under Assumptions\lll^ andl^i 

each Xq in M.'^, there exists a unique family {pt{xo,dy),t € [0,r]) of bounded 
measures on R'^ such that po{xo, .) — is the point mass at xq and 

yg G Co°°(R^R), I g{y)^{xo,dy) = J pt{xo,dy)Ltg{y). (5) 

Pt{xQ,.) is the marginal distribution at time t of the unique solution associated 
to the martingale problem for {C,C^{M.'^) ) starting from xq on [0,r]. 

Proof. Under Assumptions [l] [2] and [3] Proposition [T] implies that the martingale 
problem for £ on the domain C^{M.'^) is well-posed. For any xq in M.'^, denote 
(X, Qa;^) the unique solution of the martingale problem for C. Define 

Vte[0,r] V/eC°(R'^) Q J(xo) - E^^o [/(X*)] . (6) 

Qt is the evolution operator on [0, T] of X on C|;'(R'*). {Qt)te[o,T] is then strongly 
continuous on [0,r[. If qt{xQ,dy) denotes the law of (Xt) under Qxoj the mar- 
tingale property shows that qt{xo,dy) satisfies the equation ([5]). Integration of 
^ with respect to time yields 

qt{xo,dy)g{y) ^ g{xo) + J j qs[xo,dy)Csg{y)ds. (7) 

We have thus constructed one solution qt of ([5]) with initial condition qo{dy) = 
Cxq- This solution of ([S]) is in particular positive with mass 1. 

To show uniqueness of solutions of (O, we will rewrite ^ as the forward 
Kolmogorov equation associated with a homogeneous operator on space-time 
domain and use uniqueness results for the corresponding homogeneous equation. 
Let / S C^{R'^) and 7 G C^([0,T]) and consider the (homogeneous) dependent 
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operator A mapping functions of the form {t, x) G [0, T] x i 
will be denoted C^{[0,T]) (g) C^{R''), into : 



fi^hit), which 



(8) 



For any xq in R"^', if {X,Qxo) is a solution of the martingale problem £, then 
the law of rjt — {t,Xt) is a solution of the martingale problem for A: for any 

/eCo°°(M^O and7eC([0,T]), 



qtixo,dy)fiy)j{t)^fixohiO) 



qs{xo,dy)A{fj){s,y)ds. (9) 



Using [TTl Theorem 7.1 and Theorem 10.1, Chapter 4]), uniqueness holds for 
the martingale problem associated to the operator C on C(j"(R'^) if and only 
if uniqueness holds for the martingale problem associated to the martingale 
problem for A on CH[0, T]) (g) C^(K'^). 
Define for ah t in [0,T], for ah g in C^{[0,T]) (g, C^{R'^): 



Qtgixa) 



qt{xo,dy)g{t,y) 



(10) 



One observes that Qt-{xo) corresponds to the extension of Qt defined on the 
domain C^{R'') to the domain C'^{[0,T]) g) C^iR'^). 

Consider now a family pt{dy) of positive measures solution of ^ with po{dy) = 
Cxoidy). After integration by parts 



Pt{dy)f{y)i{t) = /(xo)7(0) + 



Ps{dy)A{f^){s,y)ds (11) 



holds true. Define for all t in [0,r], for all g in C'^{[Q,T]) ®C^{ 



oalTj>d\ 



Ptg -- 

Given ^ and (HH), for all e > 0: 

Ptih) - Peih) -- 



For any A > 0, we have 



e-^* Qt{fl){xa)dt 



Pt{dy)g{t,y). 



quixo,dy)A{fj){u,y)du 



QuiAifj))ixQ)du, 

Pu{dy)Aifj)iu,y)du^ / P^{A{fj))du. 

(12) 



fixohiO) + X 
fixohiO) + X 
fixohiO) - 



-At 



Qs{A{!-f)){xo)dsdt 



-At 



dt] Q.(A(/7))(xo)ds 



Qs{A{f^))(xo)ds 
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Consequently, 



■CXD 



e-^*Q*(A-A)(/7)(0,xo)dt = /(a;o)7(0)= / e'^' Pt{X - A){f^) dt. 



Qt defines a strongly continuous semigroup on C"([0, T] x M''). Ttie Hille-Yosida 
theorem [TTJ Proposition 2.1 and Theorem 2.6] then implies that for all A > 0, 
7^(A -A)= C^{[0, T] X M''), where 7^(A - A) denotes the image of C^{[0, T]) (g) 
C^(M'^) by the mapping g (A — A)g. Hence, since holds then for all h 
in Cj°([0,T] X M-^) 



so the Laplace transform of < H> Qth{0,xo) is uniquely determined. We will 
now show that 1 Qth (0, a:o) is right-continuous. Furthermore t — ^ b{t, .), t ^ 
a{t, .), and t — > n{t, . , .) are bounded in t on [0,T] (Assumption [T]) . It implies 
that for any fixed / G C^(R^) and any fixed 7 G Ci([0,r]), t QtA{fj){xo) 
and f — PtA{f'-f) are bounded on [0,T]. [11] Theorem 2.6] implies also that the 
set { A(/7) , / e (M'') , 7 e ( [0, T] ) } is dense in ( [0, T] x R'') and ^ shows 
that Qt-ixo) and Pf. are weakly right-continuous in t on C^{[0,T]) ®C(f (M''), 
i.e, for T >t' >t: 

hm F,,(/7) = F,(/7) hm Qt'(/7)(xo) = Qt(/7)(^o). 

Since C^([0,T] x R'') is separating [111 Proposition 4.4, Chapter 3], Qf.(xo) 
and Pj. are weakly right-continuous and (|14p holds for any A > 0, we have 
/ hiy)qt{xo,dy) = J h{y)pt{dy) for aU h e C^iR'^). This ends the proof. □ 

Remark 2.1. Assumptionsl^ [H anrfO are sufficient but not necessary for the 
well-posedness of the martingale problem. For example, the boundedness As- 
sumption [I] may be relaxed to local boundedness, using localization techniques 
developed in \2°J[ H^. Such extensions are not trivial and, in the unbounded 
case, additional conditions are needed to ensure that X does not explode (see 
JM Chapter 10]). 

2.3 Markovian projection of a semimartingale 

We will make the following assumptions, which are almost-sure analogs of As- 
sumptions [1] [2] and [31 on the local characteristics of the semimartingale ^: 

Assumption 4. 13,6 are bounded on [0,r]; 



(13) 




(14) 



3Ki > 0,yt e [0,r] ||/3t|| < Ki, \\St\\ < Ki 



a.s. 



7 



Assumption 5. /i has a density m{uj, t, dy) with respect to the Lehesgue measure 
on [0, T] which satisfies 



3X2>0,V<e[0,T] / [ihWyW) m{.,t,dy)<K2< 

JR'' 



a.s. 



and lim / m {.,t, {\\y\\ > R}) dt ^ a.s. 

Assumption 6. 

Either {i) 3e > 0, Vi G [0, T[ > e a.s. 
or (m) (5 = and there exists /3 £]0, 2[, c, i^a > 0, and a family m^(t, dy) 
of positive measures such that 
Vt G [0, T[ m{t, dy) = m'^{t, dy) + JjJ^d+p 

{lA\\yf)m''{t,dy)<K3, lim/ \\y\fm'^it,dy)=Oa.s. 

Note that Assumption [5] is only slightly stronger than stating that m is a 
Levy kernel since in that case we already have / (l A m{.,t,dy) < oo. 

Assumption [6] extends the "ellipticity" assumption to the case of pure-jump 
semimartingales and holds for a large class of semimartingales driven by stable 
or tempered stable processes. 

Theorem 2 (Markovian projection). Define, for {t, z) e [0,T] xR'^,B e BiR'^- 
{0}), 

bit,z)=E[l3t\^t- ^z], 

a{t,z)=E['dt5t\^t- ^z], (15) 
n{t,B,z) ^E[m{.,t,B)\^t- = z] . 

If (/3, (5, n) satisfies Assumptions^ O and (b, a, n) satisfies Assumption\^then 
there exists a Markov process ((Ai)tg[g , Qj^), with infinitesimal generator C 
defined by ([3|), whose marginal distributions mimick those of^: 

VtG[0,T] At =6- 
X is the weak solution of the stochastic differential equation 



Xt=^a+ [ b{u,Xu)du+ [ j:{u,Xu)dBu 
Jo Jo 

f yN{du dy)+ f f yN{du dy), 

J\\y\\<l Jo J\\y\\>l 



(16) 



where (Bf) is an n-dimensional Brownian motion, N is an integer-valued ran- 
dom measure on [0,T] x with compensator n{t,dy,Xt-)dt, N — N — n the 
associated compensated random measure and S G C°([0, T] xM.'^, Mdxn{^)) such 
that ^Y.{t,z)i:{t,z) = a{t,z). 
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We will call {X, Q_^„) the Markovian projection of ^. 
Proof. First, we observe that n is a Levy kernel : for any {t, z) e [0, T] x 



(lA n{t,dy,z)^¥. 



(lA||j/f) TO(i,dy)|6- 



< oo a.s., 



using Fubini's theorem and Assumption [S] Consider now the case of a pure 
jump semimartingale verifying (ii) and define, for B G B{M.'^ — {0}): 



m(t,dy,uj) 



cdy 



As argued above, is a Levy kernel on R''. Assumptions |4] and [5] imply that 
(6, a, n) satisfies Assumption [T] Furthermore, under assumptions either (z) or 
(ii) for ((5, m). Assumption |3] holds for {b,a,n). Together with Assumption [5] 
yields that £ is a non-degenerate operator and Proposition [T] implies that the 
martingale problem for {£t)ti£[o,T] on the domain C^(R'^) is well-posed. Denote 
{{Xt)t£[o,T], Q^o) its unique solution starting from and qti£,o,dy) the marginal 
distribution of Xt. 

Let / in C^{R'^). Ito's formula yields 

i=i Jo ClXi ^ Jo 

s<t L i=l * 

= .f{^o)+ f Vf{L-).Psds+ f Vf{Cs-).SsdWs 
Jo Jo 

+ lJ^tr[V^f{^,-)'SsSs] ds 

Vf{^,-).yMidsdy) 

(/(Cs- +y) -/(Cs-) -l{||,||<i}y.V/(e.-)) M{dsdy). 

JO JK» 

We note that 



^||y||<l 
t 



• since || V/|| is bounded /||y||<i ^fUs- )-y M{ds dy) is a square-integrable 
martingale. 

• lo I\\y\\>i '^fiL- )-y^Hdsdy) < ooa.s. since ||V/|| is bounded. 

• since V/(^5-) and Sg are uniformly bounded on [0,r], J^\/ f{£,s-)-SsdWs 
is a martingale. 
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Hence, taking expectations, we obtain: 



Jo 






ds 


s-+y)-ms-)- 


U\\y\\<i}y-^fi^s-)) M{dsdy) 




Jo 







+ K 
Observing that; 



En/(Co)]+E' 

/ / {fi^s- +y)- f{L- ) - l{|b||<i} 2/-V/(e.- )) m{s, dy) ds 
Jo Jr'^ 



<I|V/||E^ 



E^ 



< l|V^/||E 



Jo 



< DO, 



Jo 

/ / II /(^s- +y) -/(^s-) - l{||y||<i}y-V/(^s- ) II m(s,d?/)ds 

Jo Jw^ 



WSsfds 



< oo, 



/'/ 

^0 J\\v 



' m{s, dy) ds 



10 J\\y\\<l 

we may apply Fubini's theorem to obtain 



+ 2||/||E^ 



nm(s, dy) ds 
- „y||>i 



< +00, 



ds. 



= [fm + J^E^ [V/(6-)./3«] ds+\ [tr [Wfe-) ds 

+ f¥F\[ (/(e.- +y)- f{L- ) - l{||,||<i} y.V/(e.- )) m{s, dy) 

Jo iJs.'^ 

Conditioning on and using the iterated expectation property, 
= E'[/(eo)]+ rE'[V/(^,-).E'[/3.|M 

^0 

^'[tr[VVfc-)E'[*<5..5,|^,_]]]] ds 



1 

+ /*E' 



Hence 



[ ¥F [ {ms-+y)-ms-)-h\\v\\<ny-'^fi^s-)) m{s,dy)\i 
Jo I Uw^ 

= EP[/(eo)] + _^*EP[V/(6-).6(5,C.-)] ds + lj^E^ [tr[vVfe-)«(5,^.-)]] ds 
+ f I EF[{f{^s-+y)^I{^s-)'h\\v\\<i]y-yf{is-))n{s,dy,^s-)\ds. 

Jo JR'' 

/ Csf{^s-)ds 
Jo 



E'[/(6)]=E^ [f{io)]+EF 



(17) 
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Let ptidy) denote the law of (^t) under P, (|17p writes: 

/ Pt{dy)f{y)^ f Po{dy)f{y)+ f f Ps{dy)Csf{y) ds. (18) 
JW^ Jw^ Jo Jw 

Hence pt{dy) satisfies the Kolmogorov forward equation ([5]) for the operator 
C with the initial condition pa{dy) = ^o{dy) where /io denotes the law of ^o- 
Applying Theorem [H the flows gt(Co, dy) of Xt and pt{dy) of are the same on 
[0,T]. This ends the proof. 

□ 

Remark 2.2 (Mimicking conditional distributions). The construction in The- 
orem may also be carried out using 

aoit.z) =E[*5t5i|^t^ =z, Jo] , 
no{t,B,z)=W.[m{.,t,B)\it- = z, J'o] ■ 

instead of (b,a,n) in (fTSj) . If {bo,ao,no) satisfies Assumption (jS]), then fol- 
lowing the same procedure we can construct a Markov process (X, ) whose 
infinitesimal generator has coefficients (6o,ao,'^o) such that 

V/ £ C',{R%yt e [0,T] [/(6)l-^o] - E«?o [/(X,)] , 

i.e. the marginal distribution of Xt matches the conditional distribution of 
given Jo- 

Remark 2.3. For Ito processes (i.e. continuous semimartingales of the form 
([2| with fjL = 0), Gydngy U2l Theorem 4-6] gives a "mimicking theorem" under 
the non- degeneracy condition *5t.St > eld which corresponds to our Assumption 
O but without requiring the continuity condition (Assumptions^ on {b^a,n). 
Brunick & Shreve fSjj extend this result by relaxing the ellipticity condition of 
/ii^j /. In both cases, the mimicking process X is constructed as a weak solution 
to the SDE ([T6| (without the jump term), but this weak solution does not in 
general have the Markov property: indeed, it need not even be unique under the 
assumptions used in fl2^ [^. In particular, in the setting used in \12\ the law 
of X is not uniquely determined by its 'infinitesimal generator' C. This makes 
it difficult to 'compute' quantities involving X, either through simulation or by 
solving a partial differential equation. 

By contrast, under the additional continuity condition [H on the projected 
coefficients, X is a Markov process whose law is uniquely determined by its in- 
finitesimal generator C and whose marginals are the unique solution of the Kol- 
mogorov forward equation ^ . This makes it possible to compute the marginals 
of X by simulating the SDE (1161) or by solving a forward PIDE. 

It remains to be seen whether the additional Assumption\^is verified in most 
examples of interest. We will show in Section that this is indeed the case. 
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Remark 2.4 (Markovian projection of a Markov process). The term Marko- 
vian projection is justified by the following remark: if the semimartingale ^ is 
already a Markov process and satisfies the assumption of Theorem [H then the 
uniqueness in law of the solution to the martingale problem for C implies that 
the Markovian projection {X, Q^o) o/^ has the same law as (^,P^(,). So the map 
which associates (the law Q^^ of) X to may indeed be viewed as a projection; 
in particular it is involutive. 

This property contrasts with other constructions of mimicking processes fU 
\JM [TM [Sm/ which fail to be involutive. A striking example is the construction, 
by Hamza & Klebaner of discontinuous martingales whose marginals match 
those of a Gaussian Markov process. 

2.4 Forward equations for semimartingales 

Theorem [T] and Theorem [2] allow us to obtain a forward PIDE which extends the 
Kolmogorov forward equation to semimartingales which verify the Assumptions 
of Theorem [21 

Theorem 3. Let be a semimartingale given by ^ satisfying the assumptions 
of Theorem\^ Denote pt{dx) the law of onW^ . t n> pt is the unique solution, 
in the sense of distributions, of the forward equation 

Vte[0,T] ^=C:.pu (19) 

with initial condition Po = MO; where /iq denotes the law of S^q, 
where L* is the adjoint of C, defined by 

Vg G Co°°(M^R), 

Ctg{x) = -V[b{t,x)g{x)]+V^ 

+ / [g{x — z)n{t, z, X ~ z) ~ g{x)n{t, z, x) — l\\z\\<iZ.V [g{x)n{t, dz, 

where the coefficients b,a,n are defined as in (jlSp . 

Proof. The existence and uniqueness is a direct consequence of Theorem [1] and 
Theorem [21 To finish the proof, let compute £j . Viewing pt as an element of 
the dual of C^(M''), (0 rewrites : for / e C^f (M'',R) 

V/ e C„=-(M^M), J f{y)^idy) = J Ptidy)Ctf{y). 

We have 

V/ e C^{R''),Wt <t'<T < > <p,, A/ >=< £:pt,f>, 

where < ., . > is the duality product. 



a(t. x) 
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For z G R'*, define the translation operator by Tzf{x) = f{x + z). Then 
Pt{dx) Ctf{x) 



pt{dx) [b{t,x)Vf{x) + -tr [y^fix)a{t,x)] 



+ {'rzf{x)~f{x))n{t,dz,x) 

l\z\>l 



+ / {rzfix) - f{x) - z.Vfix)) nit, dz, x) 

J\z\<l 

f{x)-^[bit,x)p,idx)] + fia:)-^[^^p,{dx)] 



+ / f{x){T^z{pt{dx)n{t,dz,x)) - pt{dx)n{t,dz,x)) 

J\z\>l 

+ / f{x){T^^{pt{dx)n{t,dz,x))-pt{dx)n{t,dz,x)) 

J\z\<l 

d 

-z — {pt{dx)n{t,dz, x)) 
ox 



allowing to identify £* . 



□ 



2.5 Martingale-preserving property 

An important property of the construction of <^ in Theorem[2]is that it preserves 
the (local) martingale property: if ^ is a local martingale, so is X: 

Proposition 2 (Martingale preserving property). 



1. If ^ is a local martingale which satisfies the assumptions of Theorem\^ 
then its Markovian projection (^t)fe[o,r] local martingale on (f^o, Bt, 



2. If furthermore 



-it e [0, T] ¥F 



\\yfl^{dtdy) 



< oo, 



then {Xt)te[o,T] is a square- integrahle martingale. 

Proof. 1) If f is a local martingale then the uniqueness of its semimartingale 
decomposition entails that 



hence 



Pt+ ym{t,dy)^Q dtxF-a.e. 

A\v\\>i 

yte[Q,T], f ds b{s,Xs-)+ f yn{s,dy,X,^) 
V "'0 A\y\\>i 



= =1. 
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The assumptions on m, S then entail that X, as a sum of an Ito integral and a 
compensated Poisson integral, is a local martingale. 
2) IfE^ [J Wyfuidtdy)] < oo then 



\\yfnit,dy,Xt-) 



< oo, 



and the compensated Poisson integral in X is a square-integrable martingale. 

□ 



3 Mimicking a semimartingale driven by a Pois- 
son random measure 

The representation ([2]) is not the most commonly used in applications, where a 
process is constructed as the solution to a stochastic differential equation driven 
by a Brownian motion and a Poisson random measure 

Vie[0,r] f (3,ds+ f SsdWs+ f I ^s{y)N{dsdy), (21) 

Jo Jo Jo J 

where fo G R'', 

is a standard R"-valued Wiener process, /? and 5 are non- 
anticipative cadlag processes, is a Poisson random measure on [0,r] x 
with intensity v{dy) dt where 

/ (1 A ||2/f ) z.(dy) < oo, N = N-v{dy)dt, (22) 

and the random jump amplitude -0 : [0,r] x r2 x M'' ^ is V ® B{W^)- 
measurable, where V is the predictable cr-algebra on [0, T] x 57. In this section, 
we shall assume that 



Vte[0,r], Vf(a;,0) = and e' 



(lA||V^,(.,2/)f) v{dy)ds 

Jr<* 



< oo. 



The difference between this representation and ([2]) is the presence of a random 
jump amplitude 4't{'-^i •) in (|2T|) . The relation between these two representations 
for semimartingales has been discussed in great generality in |10[ [17] . Here we 
give a less general result which suffices for our purpose. The following result 
expresses C, in the form ([2]) suitable for applying Theorem [2j 

Lemma 1 (Absorbing the jump amplitude in the compensator). 

Ct ^ Co + f I3sds+ f SsdWs+ f I Mz)N{dsdz) 



can be also represented as 



Ct=Co+ Psds+ 6sdWs+ y M{ds dy), (23) 
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where M is an integer-valued random measure on [0,T] x R'^ with compensator 
IJ,{ui, dt, dy) given by 

VA e - {0}), n{uj, dt, A) = vii^i^ [uj. A)) dt, 

where iI}'[^{lu,A) = {z G R*^, z) G A\ denotes the inverse image of A under 

the partial map ^pf 

Proof. The result can be deduced from [lOj Theoreme 12] but we sketch here 
the proof for completeness. A Poisson random measure N on [0,T] x R"* can 
be represented as a counting measure for some random sequence (T„, C/„) with 
values in [0, T] x R'' 

N = J2hTM- (24) 

n>l 



Let M be the integer-valued random measure defined by: 



(25) 



H, the predictable compensator of M is characterized by the following property 
[161 Thm 1.8.]: for any positive P(g);B(R'')-measurable map x- [0, T] x f7 x R'* ^■ 
R+ andany AgS(R'^-{0}), 



x{s,y)M(ds dy) 



J A 



E 



X{s,y) fJ,{dsdy) 



JA 



(26) 



Similarly, for B G 6(R'^ - {0}) 



E 



X{s,y) N{dsdy) 



Uo JB 

Using formulae ([Ml) and (1^ : 



E 



x{s, y) v{dy) ds 



JB 



E 



X(s,y) M{dsdy) 



JA 



= E 



n>l 
■ t 

Jtl>7'^(..A) 
t 



10 J4>^\.,A) 

Formula (pSj) and the above equalities lead to: 

10 JA 



Y,X{TnM:Un)) 

x{s,il)s{.,z))N{ds dz) 
x{s,il)s{;z))v{dz)ds 



X{s,y)ii{ds dy) 







= E 


/'/ 

Jo Jip 



X{s,ipsi;z))v{dz)ds 



15 



Since ^/^ is a predictable random function, the uniqueness of the predictable 
compensator /i (take <j) = Id in [16, Thm 1.8.]) entails 

^j,{uj,dt,A) = v{ip^^^{iu,A))dt. (27) 

Formula (j27p defines a random measure fi which is a Levy kernel 

(lAllyf) ^i{dyds) = J^ J (lA||^,(.,2/)f) i,{dy)ds<^. 

□ 

In the case where ipt{uj, .) : invertible and differentiable, we can 

characterize the density of the compensator fi as follows: 

Lemma 2 (Differentiable case). If the Levy measure i>{dz) has a density v{z) 
and ifiptif^, .) : R'' H> R'' is a {W^ ,W^)-diffeomorphism, then C, given in ((2T]) . 
has the representation 



Ct = Q 



a + I 
Jo 



(3s ds - 



5, dW. 



yM{dsdy), 



where M is an integer-valued random measure with compensator 

m{uj;t,y)dtdy = l^^(^M''){y) IdetVyiptf^ {uj,ip^^{uj,y)) i^{ipfr^{uj,y)) dtdy, 
where V yij^t denotes the Jacobian matrix of ipt{ijJ, ■)■ 
Proof. We recall from the proof of Lemma [TJ 



E 



X{s,y)fi{ds dy) 



J A 



= E 



10 Jip7^i.,A) 

Proceeding to the change of variable ipsi-, z) = y: 

E 



x{s,'4)s{;z))v{z)ds dz 



= E 



X{s,iis{-,z))v{z)ds dz 

Ji>7\:A) 
t 

IW.CR'')} (y) y) |detVi/'sr^ (., i'^^i-, y)) y))ds dy 



J A 



The density appearing in the right hand side is predictable since is a pre- 
dictable random function. The uniqueness of the predictable compensator ^ 
yields the result. □ 

Let us combine Lemma [2] and Theorem [2l To proceed, we make a further 
assumption. 
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Assumption 7. The Levy measure v admits a density v{y) with respect to the 
Lebesgue measure on R*^ and: 

yte[0,T]3K2>0 f f {l/\Us{;y)\n Hy)dyds <K2 a.s. 

Jo J\\y\\>l 



and 



lini / 1^ i^Jt i{\\y\\ > R})) dt^O a.s. 

i?,— i-oo 



Theorem 4. Let (Q) be an Lto semimartingale defined on [0,T] by the given 
the decomposition 

Ct = Co+ f Psds+ f SsdWs+ f fMy)Nidsdy), 
Jq Jo Jo J 

where ipt[<^, .) : M'' W'- is invertible and dijferentiable with inverse (j)t{^, ■), P 
and S satisfy Assumption^ and v Assumption^^ Define 

m{t,y) = l^ye4„iR^)} IdetV^tT' (W'fe)) '^i'^^^y)). (28) 

and b{t,z) =E[l3t \Ct- = z] , 

a{t,z)=E['Stdt\Ct- ^z], (29) 

j{t, 2/, 2) = E [m{t, y) |Ct- = z] . 

If S,m satisfy Assumptions and Assumption\^ holds for {b,a,j), then the 
stochastic differential equation 

Xt^Co+ f b{u,Xu)du+ f i:{u,Xu)dB,+ [ [yj{dudy), (30) 
Jo Jo Jo J 

where (Bt) is an n-dimensional Brownian motion, J is an integer valued ran- 
dom measure on [0,T] x R** with compensator j{t,dy, Xt-) dt, J = J — j and 
E : [0,T] X M.'^ Mrfxn(IR) is a continuous function such that z)E(i, z) = 
a{t,z), admits a unique weak solution ((^t)te[o,T] i Qco) whose marginal distri- 
butions mimick those of C,: 

VtG[0,T] Xt = Ct. 
Under Q^q; X is a Markov process with infinitesimal generator C given by ([3]). 
Proof. We first use Lemma [2] to obtain the representation of 

Ct = Co+ / I3sds+ I 5sdW,+ f f yM{dsdy) 
Tfien, we observe tliat 







[ yM{dsdy) ^ ( ( yM{dsdy)+ ( ( y [M (ds dy) - fi{ds dy)] 

J Jo "'||y||<l "'0 "'||y||>l 

yM{dsdy)+ [ f yM{dsdy)- [ f y^dsdy), 

J\\y\\<l Jo J\\y\\>l Jo J\\y\\>l 
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where the terms above are weh-defined thanks to Assumption [T] Lemma [5] leads 
to: ^ ^ 

ynidsdy)^ / \\ips{.,y)\\'^ iy{y) dy ds. 

J\\v\\>l Jo "'||^a(y)||>l 

Hence: 

Ct = Co+ / Psds- [ [ Us{;y)fi^{y)dyds +f S^dWs 

Jo J\\i>s{v)\\>i J "'0 

yM{dsdy)+ f f yM{dsdy). 

J\\y\\<l Jo J\\y\\>l 

This representation has the form 1^ and Assumptions S] and [7] guarantee that 
the local characteristics of C satisfy the assumptions of Theorem [2l Applying 
Theorem [2] yields the result. □ 

4 Examples 

We now give some examples of stochastic models used in applications, where 
Markovian projections can be characterized in a more explicit manner than 
in the general results above. These examples also serve to illustrate that the 
continuity assumption (Assumption [2|) on the projected coefficients (b, a, n) in 
(|15p can be verified in many useful settings. 

4.1 Semimartingales driven by a Markov process 

In many examples in stochastic modeling, a quantity Z is expressed as a smooth 
function / : M'^ — > M of a d-dimensional Markov process Z: — f{Zt). We 
will show that in this situation our assumptions will hold as soon as Z has 
an infinitesimal generator whose coefficients satisfy Assumptions [U [2] and [3l 
Consider a time-dependent integro-difFerential operator L ~ (it)te[o.T] defined, 
for / e C^{R'^), by 

Lj(z)=5,(M).v/(z)+ x: ^^^i-) 

«j=i ' ^ (31) 

+ / [f{z + ^{t,z,y)-f{z)-ij{t,y,z).Vfiz)Hy)dy, 

jR<i 

where E : [0, T] x M'' h-^ Mdxdi^), bz : [0, T] x M'' R'' and : [0, T] x R'^ x R'' 
are measurable functions and is a Levy density. 
We assume that 

■!/'(. ,.,0) = V(^,-z, •)isaC^(R'',R'') - diffeomorphism 

3i^2>0Vte[0,r]VzeM'^ / / {l A\\^{s,z,y)\\^) v{y)dyds < K2. 

Jo "'{||y||>l} 



18 



can then be expressed as 
mz{t,y,z) = liyeM^")} IdetV?/'!"^ {t, z,il;~'^{t, z,y)) i^{ip~^{t, z,y)). 
Consider the stochastic differential equation 



Vie[o,r] Zt^Zo 



bz{u, Zu-) du 



f j:{u,Zu-)dWu 

Jo 



(32) 



'iJj{u,Zu-,y)N{du dy), 



where (Wt) is an n-dimensional Brownian motion, is a Poisson random mea- 
sure on [0, T] X M.'^ with compensator z/(y) dy dt, N the associated compensated 
random measure. Throughout this section we shaU assume that {bz, S, mz) sat- 
isfy Assumptions [U [2] and |3l Proposition [1] then imphes that for any G K'^, 
the above SDE admits a weak solution ((■Z^t)tg[o,T]: Qzo)j unique in law. Under 
Qzoj Z is a Markov process with infinitesimal generator L. Assume furthermore 
that Zt has a density qt with respect to the Lebesgue measure on R'^. 
Consider now the process 



(33) 



The aim of this section is to build in an explicit manner the Markovian Projec- 
tion of for a sufficiently large class of functions /. Before stating the main 
result, we start with an useful Lemma, which will be of importance when one 
wants to characterize the yielding Markovian projection of . 



Lemma 3. Let Z be a 

C^{R'^,R) such that 



V^-valued random variable with density q{z) and / £ 



Vz G 



(34) 



such that f{zi,--- , Zd-i, F{z)) ~ Zd- Then 
R'' M such that E[|.g(Z)|] < oo and any 



Define the junction : M — > R 
for any measurable function g : 

w e /(R^), 

E[g{Z)\f{Z)=w] 

dz,...dzd-, g(zi,---,.,-i,^)| ;<7'-------^'^<7----'--^^'-"» |. 

Proof. Consider the d-dimensional random variable k{Z), where k is defined by 

k(z) = (zi, • • • ,Zd-lJ{z)) , 

and let us compute the law of k{Z). 

/ 1 \ 

^ ' ••• 1 

, M. ... df df_ I 

\ dzi dzd-i dzd / 
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One observes that |det(V2K;)|(2;) = §^{z) 
diffeomorphism with inverse k~^. 



> 0. Hence k is a 



k{k ^{z)) = (k/(2;),--- ,K^^^{z),f{Ki\z),--- ,k/{z)) = z. 

For \ <i< d— 1, = Zi and /(-Zi, • ■ • , -Zd-i, k;^^(^;)) = z^ that is ^^^(-z) = 

Hence 

Define ^^(-z) c^-z the inverse image of the measure q{z) dz under the partial map 
K by 

Qk{z) = l{^(jsid)}iz)\detC^,K~^)\{z)q{K~'^{z)) 



df 



dzd 



{zi,-- - , Zd-i,F{z)) q{zi,- ■ ■ , Zd-i, F{z)). 



qniz) is the density of k{Z). So, for any w G 
E[g{Z)\f{Z)=w] = E[g{Z)\K{Z) = {zu--- ,Zd-i,w)] 



/ 



dzi...dzd-ig{zi, - ■ ■ ,zd-i,w) 



q{zi,--- ,Zd-i,F{zi,--- ,Zd-i,w)) 



Of 



dzd 



(2:1, • • ■ , Zd-1, F{z-i_, ■■■ , Zd-1, w)) 



□ 



We can now formulate the main result of this section: 
Theorem 5. Let f G C^(M'^,M) with bounded derivatives such that 



yz G M<* 1^(2) ^ 0. 

OZd 

Define, for w e /(M''),t e [0,T], 

b{t, w)= f \vf{.).hz{t, .) + itr [VV(-)*S(i, .)] 

+ / {f{.+xl;{t,.,v))- f{.)-xl;{t,.,y).Vf{.)) v{y)dyMzu--- ,Zd-uw) 



(35) 



K<1 



<j{t, w) 



(zi,-- - ,Zd-i,F{zi,--- ,Zd-i,w)) 
[ ||V/(.)S(i,.)ll'(^i,--- ,Zd-i,w) 

jRd-l 



qt{zi,--- ,Zd-i,F{zi,--- ,Zd-i,w)) 



df. 

dzd 



{zi,--- ,Zd-i,F{zi,--- ,Zd-i,w)) 



1/2 



(36) 
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and the measure j{t, du, w) defined on M — {0} by 
j{t,[u,oo[,w) = i l{f^,+^(t^.^y))-f(.)>u}{zi,--- ,Zd-i,w)i'{y)dy 



,Zd-i,F(zi, ■ ■ ■ ,Zd^i,w)) 



§^{zi,--- ,Zd-i,F{zi,--- ,Zd-i,w)) 



(37) 



for u> and analogously for u <0. Then the stochastic differential equation 



— ^0 



b{s,X^)ds+ I a{s,X,)dBs 

/ / yJ{dsdy), 

'o A\y\\>i 



(38) 



yj{dsdy) 

'0 J\\y\\<l 

where (Bt) is a Brownian motion, J is an integer-valued random measure on 
[0, T] X K with compensator j{t, du, Xt-)dt, J = J — j, admits a weak solution 
{{^t)te[OT],Qso)' unique in law, whose marginal distributions mimick those of 

vte[o,T] Xt=^f 

Under Q^^, X is a Markov process with infinitesimal generator £ given by 



V/ e Co°°(R),A/(w) = b{t,w)f'{w) + 



■f"H 



+ [ [f{w + u)- f{w)-uf'{w)]3{t,du,w). 
Proof. Applying Ito's formula to f{Zt) yields 

6 



= 6+/ Vf{Z,-).bz{s,Zs-)ds+ [ Vf{Z,-).i:{s,Z,_)dWs 
Jo Jo 

+ 11 tr[VV(^.-)*S(s,Z,_)E(s,Z,_)] ds + J^ Vf{Z,-)4{s,Z,_,y) N{dsdy) 
+ f I {f{Z,-+i;is,Z,_,y))-f{Z,-)-i;{s,Z,_,y).Vf{Z,-))N{dsdy) 

Jo JRt' 

= ^0 + ^ [v/(Z,-).6z(s,Z,_) + itr[vV(^«-)*S(s,^,_)E(s,Z,_)] 

+ / {f{Z,-+i;is,Z,-,y))-f{Z,-)-ij{s,Zs-,y).WfiZ,-))u{y)dy 

J»A 

+ f\f{Z,-).i:{s,Z,_)dWs + f [ {f{Z,- +ilj{s,Z,_,y))-f{Z,-)) N{dsdy). 

Jo Jo JM<* 



If S satisfies assumption {i) then (-Bt)tg[o,T] defined by 

_ vfjz, ).s(/.z,_)n-, 

* ||V/(Z,_)E(i,Z,_)|| ' 
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is a continuous local martingale with [B]t = t. Levy's theorem implies that B 
is a Brownian motion. 

If E = and (ii) holds, then is a pure-jump semimartingale. Define K-t 

with ^(i, z, y) — ip{t, z, Hz{y)) where 
K^(y):M'^ ^ 

y ^ ivi,--- ,yd-^iJ{z + y)- f{z)). 

Since for any z G Mf^, |detVj,Kz| (y) = 

i^:\y) = {yir-- .yd-i,F4y)) F,iy):R''^R • • • , y^.i, F,(y)))-/(z) = y^. 

Considering the inverse of "0 that is ip{t, z, y), z) — y, define 

<i>it,z,y)=<j){t,z,K,-\y)). 

$ corresponds to the inverse of 4* and $ is differentiable on K'' with image M''. 
Define the random measure m(t, z, y) dt dy 

m{t,z,y) = \delVy^{t,z,y)\ v{^{t,z,y)) 
= \<leiVy4){t,z,K';^{y)) 
Using Assumption [71 

/ / {l^\\^{s,z,y)f) v{y)dyds 

Jo J\\y\\>l 

- r / (lA(V/(s,z,y)2 + ... + V'''-i(s,z,2/)2 + (/(z + ^(s,z,2/))-/(z))2)) j.(y)dyds 
Jo J\\y\\>l 

< I I {lh{^\s,z,yf + --- + i;''-\s,z,yf + \\VffU{s,z,y)f))v{y)dyds 
Jo J\\y\\>l 

< f I {I ^ {2y \\V fr)U(s, z,y)f) v{y)dyds 
Jo J\\y\\>l 

is bounded. One may apply Lemma [5] and express ICt as K-t ~ /J J yM{ds dy) 
where AI is a compensated integer-valued random measure on [0,r] X M'' with 
compensator m(t, Zt-, y) dy dt. 

Extracting the d-th component of /Ct, one obtains the semimartingale decom- 
position of on [0, T] 

it=io+ I Psds+ I 5sdB,+ [ f uk{dsdu), 
Jo Jo Jo J 



a 

8yd 



fiz + y) 



> 0, one can define 



§^^{z + .-Hy)) 



v{(j){t,z,K^^{y))). 
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where 

'Pt - Vf{Z^-).bz{t, Zt-) + itr [v2/(^t-)*S(t, Zt_)S(t, Zt^)] 

+ 4, {f{Z,- + ^{t, Zt-,y)) - fiZ,- ) - ijit, Zt-,y).Vf{Z^- )) u{y) dy, 
6t ^\\VfiZt-Mt,Zt-)\\, 

and K is an integer- valued random measure on [0,r] x M with compensator 
k{t, Zt-,u) du dt defined for ah z G M'' via 



k{t,z,u)^ / m{t,{yi,--- ,ya-i,u),z)dyi--- dyd-i 

= / \deWy^{t,z,{yi, - ■ ■ z, (yi, • ' • ,yd-i,u))) dyi ■ ■ ■ dyd-i, 

(39) 

and K its compensated random measure. In particular for any u > 0, 

k{t, z,[u, oo[) = l{f(z+4,(t,z,y))^f(z)>u}v{y)dy. (40) 



Let us show that if (6z,I],i/) satisfy Assumptions [l] [2] and |3] then the triplet 
{6t, Pt,k{t, Zt-,u)) satisfies the assumptions of Theorem[2] First, note that /?* 
and 6t satisfy Assumption |4] since bz{t,z) and I]{t,z) satisfy Assumption [T] and 
V/ and V^/ are bounded. 



{lA\\u\\^) k{s,Zs-,u)duds = I \ {lA\f{Z,- +^s,Zs-,y))- f{Z,-)\^) v{y)dyds 

(lA||V/||2||^(s,Z,_,w)f) v{y)dyds, 



< 



is bounded by Assumption [T) Hence k satisfies Assumption [5] 

As argued before, one sees that if S is non-degenerate then St is. In the case 

St = 0, for t G [0,T[,R > 0,z G B{0,R) and g G Cg(R) > 0, consider C and 
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Kt > chosen via Assumption [3] 

k{t, z, u) 



\deWy^{t,z,(yi, - ■ ■ ,yd-i,u))\ z, (yi, • • • ,yd-i,u))) dyi ■■■ dyd^ 



\detVy(l>{t, z, ^(yi, • • • , y^^i, 



^{z + K^^{yi,--- ,yd-i,u)) 
dyd 



v{(j){t, z, • • • , yd, u))) dyi ■■■ dyd- 



> 



dyd 



{z + K^^{yi,--- ,yd-i,u)) 



C 



c 



-1 ,2/d-i,«)r+^ 

1 f c 



ll't^ (yi, • • • ,yd-i,u) 

dyi ■■■ dyd-i 



■dyi ■ ■■ dyd- 



1 



\\{y^/u,--- ,yd-i/u,l)r+P 



^dyi--- dyd-i 



C 



u 



|l+/3' 



with C" = J^a-i C \\{wi, ■ ■ ■ ,Wd-i, 1)11 ^ dwi ■ ■ ■ dwd 
Similarly 



(1 A [k{t,z,u) 



c 



du 



(lA|u|^) 



dl 

dyd 



(z + K^^(yi,--- ,yd-i,u)) 



|detVy(/)(t, z, ,yd-i,u))\ t/((/)(t, z, ^(yi, ■ • • ,yd^i,u))) 



C 



dyi ■■■ dyd-i du 



< 



(lA|/(z + (yi,-.- ,yd_i,«))-/(z)|'5) 

(^\detVy(t){t,z,{yi,--- ,yd-i,u))\ v{4>{t,z,{yi,- ■ ■ ,yd-i,u))) 

~TT7 ~ \iidTd) ■ • • dyd-i du 

11(2/1, • • • ,2/d-i,u)r+'^/ 

(1 A ||V/||I|(yi,--- ,yd-i,u)f) ( \deWy(t>{t, z, {yi, ■ ■ ■ ,yd-i,u))\ iym,z,{y 



C 



dyi ■■■ dyd-i du 
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is also bounded. Similar arguments would show that 

C 



lim 



«|<e 
T 



\u\P [k{t,Zt-,u) 



|l+/3 



du = a.s. 



and lim / k [t, Zt-,{\u\> R}) dt ^ Q a.s., 

R->-oo Jo 

since this essentially hinges on the fact that / has bounded derivatives. Define 
as in Theorem [5] 

6(i,u.) -E[A|/(Zt_) = 
ait,w)=E[Sf\fiZt-) = w]'^\ 
jit,u,w) =E[kit, Zt-,u)\f{Zt-) ^ w] . 

Applying Lemma [3l one can compute explicitly the conditional expectations 
above. For example, 



b{t,w) 



V/(.).&z(t,.) + -tr [v2/(.)'S(i,.)S(t, 



+ / (/(. + y)) - /(■) - ^(i, 2/)-V/(.)) i^iy) dy, (zi, 
qt{zi,--- ,Zd-i,F{zi, - ■ ■ ,Zd-i,w)) 



.Zd-i,F{zi, - ■ ■ ,Zd-i,w))\ 



with F : R'^ ^ R defined by f{zi, ■ ■ ■ ,Zd^i,F{z)) z^. Since / is with 
bounded derivatives, {bz,^,v) satisfy Assumption [2] and {t,z) qt{z-) is con- 
tinuous in (i, z) on [0,T] x R'^ (see Theorem [T]) then h{t, .) is continuous on R 
uniformly in t E [0,T]. Hence Assumption [2] holds for b. Proceeding similarly, 
one can show that that Assumption [2] holds for a and j so Theorem [2] may be 
applied to yield the result. □ 



,Zd-l,w) 



4.2 Time changed Levy processes 

Models based on time-changed Levy processes have been the focus of much 
recent work especially in mathematical finance [6]. Let Lt be a Levy process, 
(6, cr^,i^) be its characteristic triplet, N its jump measure. Define 

Jo 

where (6*4) is a locally bounded J^t-adapted positive cadlag process, interpreted 
as the rate of time change. 

Theorem 6 (Markovian projection of time-changed Levy processes) . Let L be 
a scalar Levy process with triplet (6, cr^,i^) and let ^4 = Lq^ with Qt — j^^sds 
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where 9t > is a positive semimartingale. 
Define 

yt G [0, T] \/zgR a{t, z) = E[6t\it- = z], 

and suppose that a{t,.) is continuous on R'^, uniformly in t G [0, T]. Assume 
that 

lim > i?}) = 

and either {i) or {ii) holds for {cr,9t y), then 

• (Ct) has the same marginals as {Xt) on [0, T], defined as the weak solution 
of 

Xt=^o+ I (J\/a{s,Xs-)dBa + [ ba{s,Xs-)ds 
Jo Jo 

+ / zJ{dsdz)+ / / zJ{dsdz), 
Jo J\z\<i Jo J\z\>l 

where Bt is a real-valued brownian motion, J is an integer-valued random 
measure on [0,T] x R with compensator a{t,Xt-) y{dy) dt. 

• The marginal distribution pt of is the unique solution of the forward 
equation: 

dpt 
dt 

where, is given by 



Ct-Pu 



C*tg{x) = -b— [a(t, x)g{x)] + Y^[a'(<, x)g{x)] 



+ 



[ v{dz) 
Jw 



g{x - z)a{t,x- z) - g{x)a{t,x) - l\\^\\<iz. — [g{x)a{t,x)] 



with the given initial condition po{dy) = jJ-oidy) where /zq denotes the law 

Proof. Consider the Levy-Ito decomposition of L: 

Lt=bt-haWt+ I I zN{dsdz)+ I / zN{dsdz). 

Jo J\z\<l Jo J\z\>l 



Then ^ rewrites 



Ct = Co + (yW{&t) + bOt 

ret r rOt 



+ f ' f zN{dsdz)+ I ' j zN{dsdz). 

Jo J\z\<l Jo J\z\>l 
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W(9f) is a continuous martingale starting from 0, with quadratic variation 
Qt = JoSgds. The Dubins-Schwarz theorem (see [HI Theorem 1.7]) implies 
that one can pick Z a Brownian motion independent of W, such that 

w{et)= f x/FsdZs. 

Jo 

Hence Xt is the weak solution of : 

Xt = Xq+ [ ayfd.dZ, + / hds ds 
Jo Jo 

+ [ f z9sN{dsdz)+ [ [ zdsN{dsdz). 

Jo "'|z|<l Jo "'Iz^l 

Using the notations of Theorem [2l 

Pt^bOt, St^cjy/o't, m{t,dy) ^ 6tv{dy). 
Since Assumptions |4] and [5] are satisfied and : 

b{t,.)=¥.[f5t%- ^ ]^ba{t,.), 

n{t,B,.)=¥.[m{t,B)\i^- = ] ^ a{t, :)v{B), 

are all continuous on R uniformly in t on [0, T]. One may apply Theorem [21 □ 

The impact of the random time change on the marginals can be captured 
by making the characteristics state dependent 

( ba{t,Xt-),a'^a{t,Xt-),a{t,Xt^)u ) 

by introducing the same adjustment factor a{t^Xt-) to the drift, diffusion co- 
efficient and Levy measure. In particular if a{t, x) is affine in x we get an affine 
process [8] where the affine dependence of the characteristics with respect to the 
state are restricted to be colinear, which is rather restrictive. This remark shows 
that time-changed Levy processes, which in principle allow for a wide variety of 
choices for 9 and L, may not be as flexible as apparently simpler affine models 
when it comes to reproducing marginal distributions. 
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